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Abstract 

We perform a new calculation of the hadronic contributions, a(Hadronic) 
to the anomalous magnetic moment of the muon, a^. For the low energy contri- 
butions of order we carry over an analysis of the pion form factor F,r (t) using 
recent data both on e'*'e~ tt^tt" and VrT^^T^'^ ■ In this analysis we take 

into account that the phase of the form factor is equal to that of tttt scattering. This 
allows us to profit fully from analyticity properties so we can use also experimental 
information on at spacelike t. At higher energy we use QCD to supplement 

experimental data, including the recent measurements of e^e~ hadrons both 
around 1 GeV and near the cc threshold. This yields a precise determination of 
the O(a^) and 0(0?) -I- 0(o?') hadronic part of the photon vacuum polarization 
pieces, 

lO" X oF' (h.v.p.) = 6 909 ± 64; lO" x a'^+^' (h.v.p.) = 7 002 ± 66 

As byproducts we also get the masses and widths of the p", and very accurate 
values for the charge radius and second coefficient of the pion. 

Adding the remaining order o? hadronic contributions we find 

lO" X a*°°'^ (Hadronic) = 6 993 ± 69 (e+e" + r + spaccl.) 

The error above includes statistical, systematic and estimated theoretical errors. 
The figures given are obtained including r decay data; if we restrict ourselves to 
e~ data, slightly lower values and somewhat higher errors are found. 

This is to be compared with the figure obtained by subtracting pure 
electroweak contributions from the recent experimental value, obtained from mea- 
surements of the muon gyromagnetic ratio — 2), which reads 

10" X a'^'^P- (Hadronic) = 7 174 ± 150. 
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1. Introduction 

The appearance of a new, very precise measurement of the muon magnetic moment'^l has triggered the 
interest in theoretical calculations of this quantity. Particularly, because the experimental figure (we give 
the result for the anomaly, averaged with older determinations!^!) 

10" X a^(exp.) = 116 592 030 ± 150 (1.1) 

lies slightly above theoretical evaluations based on the standard model, as much as 2.6cr in some cases. 

It should be noted that all modern^ theoretical determinations''^"^] are compatible among them- 
selves within errors (of order 100 x 10~^^) and that, with few exceptions, they are also compatible with 
the experimental result, (1.1), at the level of 1.5cr or less. Because of this, it is our feeling that a new, 
complete evaluation would be welcome since, in fact, there exists as yet no calculation that takes fully into 
account all theoretical constraints and all the new experimental data. These experimental data allow an 
improved evaluation of the low energy hadronic contributions to a^, both directly from e+e~ annihilations 
(in the p regionl^l and around the </> resonance'^!) and, indirectly, from r decays'^^' and, also indirectly, 
from measurements of the pion form factor in the spacelike region. l^^l Moreover, the BES'^^l data, covering 
e+e" annihilations in the vicinity of cc threshold, permit a reliable evaluation of the corresponding hadronic 
pieces. In fact, the main improvements of the present paper arc the calculation of the two pion contribution 
to the hadronic part of a^^, using all available experimental information and fulfilling compatibility with all 
our theoretical knowledge, and the pinning down of the multipion, KK and cc contributions. This we do in 
Sects. 3, 4 (in Sect. 2 we formulate the problem). In Sect. 5 we discuss other hadronic corrections, including 
one that, as far as we know, has been hitherto neglected, and which, though small (~ 46 x 10~^^) is relevant 
at the level of accuracy for which we are striving. Finally, in Sect. 6 we discuss our results and compare 
them with experiment. 

The main outcome of our analysis is an accurate and reliable determination of the hadronic contribu- 
tions to at order a^. In fact, in all regions where there are difficulties we perform at least two evaluations, 
and take into account their consistency (or lack thereof). Furthermore, we discuss in some detail (including 
ambiguities) the 0{a^) hadronic contributions. 

As a byproduct of the low energy calculations we can also give precise values for the p+ masses 
and widths. 



mpo = 772.6 ± 0.5 McV, Fpo = 147.4 ± 0.8 MeV; 
mp+ = 773.8 ± 0.6 MeV, = 147.3 ± 0.9 MeV; 



(1.2) 



(1.4) 



for the P-wave tttt scattering length, 

al = (41 ± 2) X 10"^ m"^, (1.3) 

and for the pion mean squared charge radius and second coefficient: 

{rl) = 0.435 ± 0.002 fm^ c„ = 3.60 ± 0.03 GeV"^ (e+e" + t + spacelike). 

{rl) = 0.433 ± 0.002 fm^, c„ = 3.58 ± 0.04 GeV"^ (e+e" -|- spacelike). 

We give results both using only direct results on F„, from e+e" annihilations, or involving also the decay 
r+ Pi-tt+tt'', which last we consider to be our best estimates. 
So we write, 

mil m J • \ f 6 993 ±69 (e+e~ + t + spacel.) ,^ 
10 xa(Hadromc) = |gg^3^gg (e+e" + spacel.) ^^'^^ 

Note that in a(Hadronic) we include all hadronic contributions, 0{a^) as well as 0{a^). The errors 
include the statistical errors, as well as the estimated systematic and theoretical ones. This is to be compared 
with the value deduced from (1.1) and electroweak corrections 

10" X o'^'^P- (Hadronic) = 7174± 150, 

from which (1.5) differs by l.lcr. 



^ By modern we here mean, somewhat arbitrarily, those obtained since 1985. A more complete list of references, 
including earlier work, may be found in ref. 7. 
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2. Contributions to 

We divide the various contributions to as follows: 

= a(QED) + a(Weak) + a(Hadronic). (2.1) 

Here a(QED) denote the pure quantum electrodynamics corrections, and a(Weak) are the ones due to W, Z 
and Higgs exchange. The hadronic contributions can, in turn, be split as 

a(Hadronic) = a^'^\h.Y.p.) + a(Other hadronic, 0{a^)). (2.2) 

a*-^-* (h.v.p.) arc the corrections chie to the hadronic photon vacuum polarization contributions (Fig. 1), 
nominally of order (see Subsects. 3.3 and 5.2 for a qualification of this statement). We will discuss in 
detail the "Other hadronic, 0(a^)" in Sect. 5. 

According to the review of Hughes and Kinoshitat^^l one has 

10^^ X a(QED) = 116 584 705 ± 1.8 

(2 3) 

10^^ X a(Weak) = 151 ± 4. 



There is no dispute about these numbers. If we combine them with (1.1), we can convert this into 
a measurement of the hadronic part of the anomaly: 

10" X a'^^P- (Hadronic) = 7 174 ± 150. (2.4) 

Our task in the present paper is the evaluation of this quantity. 

We now say a few words about the piece a^^^(h.v.p.), which is the most important component of 
a(Hadronic). As Brodsky and de Rafael^^^l have shown, it can be written as 



/>oo 

a(2)(h.v.p.) =127r / dtK{t)Imn{t) 

J Ami 



4m?. 



Kit) =5-5TK(i); K{t) = / dx 



-L 



"inH ^ " ^ ' Jo a;2 + (1 - x)t/ml' 

Here U is the hadronic part of the photon vacuum polarization function. An alternate formula is obtained 
by expressing Im U in terms of the ratio 

/•CO 

a(2)(h.v.p.) = / dtK{t)R{t). (2.4b) 
The superindex (0) here means 'lowest order in the electromagnetic interactions'. 
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At low energy {t < 0.8 GeV^) we can separate the contribution from three pion states and that from 
two pions. The first will be discussed in Sect. 4. The two pion contribution in turn can be expressed in 
terms of the pion form factor, F^^, 

Imi72.(t) = — \F^{t)\\ (2.5) 

so that, for the two-pion contribution up to energy squared t^, 

fto / 4^2 \ 3/2 

a^{2n;to) = \j^ " mW^t)?- (2.6) 



3. The pion form factor 
3.1. Theory 

The evaluation of the pion form factor is slightly complicated by the phenomenon of a; — p interference. This 

can be solved by considering only the isospin 1=1 component, and adding later the ^ 27r and interference 
in the standard Gounaris-Sakurai way. This is equivalent to neglecting, in a first approximation, the breaking 
of isospin invariance. We will also neglect for now electromagnetic corrections. In this approximation the 

properties of F-,^{t) arc the following: 

(i) FTjit) is an analytic function of t, with a cut from 4to^ to infinity. 

(ii) On the cut, the phase of F.^{t) is, because of unitarity, identical to that of the P-wave, 7 = 1, tttt 
scattering, 5\{f?}, and this equality holds until the opening of the inelastic threshold ai t = to (Fermi- 
Watson final state interaction theorem). 

(iii) For large t, F„{t) ~ 1/f. Actually, one knows the coefficient of this behaviour, but we will not need it 
here. 

(iv) F(0) = 1. 

The inelastic threshold occurs, rigorously speaking, at t = 16m^. However, it is an experimental 
fact that inelasticity is negligible until the quasi-two body channels coir, anr . . . are open. In practice we will 
take 

to ~ 1 GeV^ 

and fix the best value for to empirically. It will be to = 1-1 GeV^, and we will see that, if we keep close to 
this value, the dependence on to is very slight. 

The properties (i-iv) can be taken into account with the well-known Omnes-Muskhelishvili method. 
We construct a function J{t) with the proper phase by defining 

Jit) =ex.p{- ds \ +- ds \ } . (3.1a) 

We have written the dispersion relation with one subtraction to ensure that J(0) = 1. The singular integrals 
are understood to be calculated replacing t ^ t + ie, e > 0, and letting then e ^ 0. In particular, we have 

|J(i)|=exp(*P.P. r ds-^ + ^ Tds-^l, Aml<t<to. (3.1b) 

Defining then the function G by 

F^it) = G{t)J{t), (3.2) 

it follows from properties i-ii that G{t) is analytic with only the exception of a cut from to to infinity, as we 
have already extracted the correct phase below t = to- 

We can, in Eq. (3.1a), take any value we like for the phase Sl{t), as a change of it only results in 
a redefinition of G; but it is convenient to choose 6l{t) so that it joins smoothly Sl{t) at t = to to avoid 
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spurious singularities that would deteriorate the convergence, and so that J has the correct behaviour at 
infinity. Both properties are ensured if we take, simply, 

so that Sl{to) = Sl{to) and, for large t, 5l{t) w and we recover the behaviour 1/t of the form factor. Then 
we can rewrite more explicitly (3.1) by integrating the piece with 

(,_i) e.p{i/^^a,Al£Lj. ,,3, 

3.1.1. The phase si 

We can apply the effective range theory to the phase According to this, the function 

^it)^^^cotSl{t), fc=^/*HM (3.4a) 

is analytic in the variable t except for two cuts: a cut from — oo to 0, and a cut from t = to to +oo. To profit 
from the analyticity properties of tjj we will make a conformal transformation.^ We define 

^^ Vt-Vh^ _ (3.4b) 

When t runs the cuts, w goes around the unit circle. We may therefore expand tp ina. power series convergent 
inside the unit disc. However, the existence of the p resonance implies that we must have coti5j(TOp) = 0. 
It is therefore convenient to incorporate this piece of knowledge and expand not tjj itself but the ratio 

'ij){t)/{TV?p — t)= ip{t): so we write, 

tl^{t) ^ {ml - t)t}){t) = {ml - t) {bo + b^w + 62^^ + •••}. (3.4c) 
The P-wave, / = 1 tttt scattering length,"^ a\, is related to i/' by 

= mJiAmir (^-^^^ 

Likewise, from the relation 

1 const. 



^p = -7T7^' K = - Ami. (3.5b) 



cot ,51 {t) - i t^ml ml-t- 2m/t^/^{t) 
we find the expression for the rho width: 

TOpV'(TOp) 

Experimentally, [^''l a\ ~ (0.038 ± 0.003)m~^, and, according to the Particle Data Group Tables, t^®' mp = 
769.3 ± 0.8; Fp = 150.2 ± 0.8 MeV. Note, however, that we do not assume the values of Wp, Fp. We only 
require that ij) has a zero, and will let the fits fix its location and residue. 

It turns out that, to reproduce the width and scattering length, and to fit the pion form factor as 
well (see below) , only two terms in the expansion are needed, so we approximate 

*i (,) = Arc CO. { - + --l^^] } ; (3^6) 



^ The method of conformal transformations is rigorous, simpler and produces better results than that employed in 
ref. 4. 

^ For details on tttt scattering, including analyticity properties and the Fermi-Watson theorem, see e.g. ref. 15. More 
details on the solution of the Omnes-Muskhelishvili equation can be found in N. I. Muskhelishvili, Singular Integral 
Equations, Nordhoof, 1958. 



-4- 



-PRECISION DETERMINATION OF THE PION FORM FACTOR AND CALCULATION OF THE MUON Q — 2- 



TUp, bo, bi are free parameters in our fits. 
3.1.2. The function G{t) 

Because we have already extracted the correct phase up to t = to, it follows that the function G{t) 
is analytic except for a cut from t = to to +00. The conformal transformation 



(3.7a) 



5\/io + Vto — t 

maps this cut plane into the unit circle. So we may write the expansion, 

G{t) = l + Ao + ciz + C2Z^ + C3Z^ + ■■■ (3.7b) 

that will be convergent for all t inside the cut plane. We can implement the condition G(0) = 1, necessary 
to ensure -F,r(0) = 1 to each order, by writing 

^0 = - [ciZo + C2Zq + C3Z0 H ] , zo = z\t=o = -1/3. 

The expansion then reads. 



Git) = l + ci{z + 1/3) + C2(z2 - 1/9) + cs{z^ + 1/27) + • • • . (3.8) 
We will need two-three terms in the expansion, so we will approximate 

Ci, C2 free parameters. 

An interesting feature of our method is that, even if we only kept one term in each of the expansions 

(3.6, 8), that is to say, if we took &i = ci = C2 = 0, we could reproduce the experimental data with only a 
15% error; so we expect (and this is the case) fast convergence of the series. It is important also that our 
expression for F.n-(f) is valid in the spacelike as well as in the timelike region, provided only t <to- What is 
more, (3.6, 8) represent the more general expressions compatible with analyticity, the Fermi-Watson theorem 
and the effective range theory, which follow only from the requirements of unitarity and causality. Therefore, 
by employing our expansions, we do not introduce uncontrolled biases in the analysis, and hence we minimize 
the model dependent errors.^ 



G{t) = 1 + ci 



to + \/to — t 



3.2. Fits 

In order to fit Fjr, and hence get the 27r low energy (4m^ <t< 0.8 GeV^) contribution to a^^^(h.v.p.), we 
have available three sets of data: 

• e+e~ 7r+7r~, t timelike (Novosibirsk, ref. 8). 

• F^{t), t spacclike (NA7, ref. 11). 

• In addition, one can use data from the decay r+ Vt-tt^tt^ (Aleph and Opal, ref. 10). 

For this last we have to assume isospin invariance, and neglect the isospin 7 = 2 component of tt+tt^, 
to write the form factor v\ for r decay in terms of F^^: 

^1 = 1^(1-^) ' l^-WI'. (3-9a) 
where, in terms of the weak vector current = u'^^d, and in the exact isospin approximation, 

= {-p'g^^.+Pi.Pu) n^{t) = l / d^xe">-^0\TV+{x)V,{0)\0); VI = 27rlm77^. (3.9b) 



* The remaining approximations are neglect of the inelasticity between I6rn^ and to, experimentally known to be at 
the 10~^ level or below, and we have the errors due to the truncation of the expansions; we will also check that 
they are small. By contrast, other functional forms used in in the literature are either incompatible with the phase 
of Fjr, or with its analyticity properties (or both), which causes biases in the fits. 
The errors due to breaking of isospin and electromagnetic corrections will be discussed below. 
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Before presenting the results of the fits a few matters have to be discussed. A first point to clarify 
is that we will not include in the fits the old data on F^r in the spacelike or timelike regions, or on pion-pion 
phase shifts'^^l. We have checked that, if we add the first two sets, the results of the fit vary very little (see 
below); but they cause a bias. This is so because there is incompatibility^ between old spacelike and timelike 
data, and also with data on nn phase shifts, already noticed in CLY.!"*! Doubtlessly, this is due to the fact 
that most old data for spacelike momentum were extracted from processes with one pion off its mass shell 
so that models were necessary to extrapolate to the physical form factor. In fact, a very important feature 
of the NAT^^^l data is that they are obtained from scattering of real pions off electrons, hence we do not 
require models to extract F^^ from data. 

The reason for the model dependence of tttt phase shift analyses is that these are extracted from fits 
to ttN —> ttttN scattering and thus require a model for the pseudoscalar form factor of the nucleon, a model 
for the interactions of the nucleon and the final state pions, and a model for the dependence of tttt scattering 
on the mass of an external pion. Indeed, different methods of extrapolation result in different sets of phase 
shifts, as can be seen in the experimental papers of Hyams et al. and Protopopescu et al., ref 17, where five 
different determinations are given. However, we will use the scattering length, a\, and employ the tttt phase 
shifts as a very important a posteriori test of our results. 

We could consider, besides this information, to include as input the values of several quantities that 
can be estimated with chiral perturbation theory methods, such as (r^) and a\. We do not do so because 
the problem with these calculations is the estimate of their errors, a difficult matter; so we have preferred 
to avoid possible biases and instead obtain these quantities as byproduct of our calculations. Then we check 
that the results we get for all of them are in agreement, within errors, with the chiral perturbation theory 
results; see below. With respect to a\ we actually constrain it to the region obtained from tttt scattering 
experimental data only; its error is chosen such that it encompasses the various values given in the different 
experimental determinations (ref. 17). We take, 

al = (38 ±3) X 10-3 m-3; 

we will see that the value our best fit returns for this quantity is satisfactorily close to this, as indeed we get 
(41 ±2) X 10-3 m-3. 

Another remark concerns the matter of isospin breaking, due to electromagnetic interactions or the 
mass difference between u, d quarks, that would spoil the equality (3.9a). It is not easy to estimate this. 
A large part of the breaking, the lo —^ 2t: contribution and uj — p mixing, are taken into account by hand, 
but this does not exhaust the effects. For example, merely changing the quark masses from + ni-,^- to 
TO^o + m^o, in a Breit-Wigner model for the p, shifts a(^^(h.v.p.) by ~ 50 x IQ-"'^'^, so a deviation of this 
order should not be surprising.^ 

As stated above, Eqs. (3.9) were obtained neglecting the mass difference to„— mrf and electromagnetic 
corrections, in particular the tt*^ — tt"*" mass difference. We can take the last partially into account by 
distinguishing between the pion masses in the phase space factor in (3.9a). To do so, write now (3.9b) as 

77^; = i j d''xe'P-^Q\TV+{x)V,{Qm = {-p" 9t.. + Pf.P.) [t) + p^p.n^ ; vi = 27rlm77^. (3.10a) 
We get 

rr im^,-m^.)^]\^ ("^'^-+--'')'n'^'|7^,(t)P. (3.10b) 



^1 = 1^ 



t 



t 



To compare with the experimentally measured quantity, which involves all of Im 7T^^ , we have to neglect the 
scalar component 77'^, which is proportional to (m^ — m„)^, and thus very small. 

To understand the situation we will proceed by steps. First of all, we start by fitting separately 
e+e^ and t data, in the timelike region, using (3.9a) (we remark that although in a(27r; t < 0.8 GeV^) only 



^ At the level of 1.5 to 2a. 

^ The relevance o f isospin breaking in this context was pointed out by V. Cirigliano, G. Ecker and H. Neufeld, 



hep-ph/0104267, 2001 
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enter the values of FttC^) for 4mJ to 0.8 GeV^, we fit the whole range up to f = = 1.1 GeV^). Then, we 
get quite different results: 

.^n8nv2^ / 4 715 ±67 (e+e"; xVd-O.f. = 106/109 = 0.96) , , . 

a(2,r;t<0.8GeV) = |^g^^^26 (r; X^d-cf. = 52/48 i 1.09). ^^'^^"^^ 

This takes into account statistical errors only for e+e", but includes systematic ones for r decay as these 
are incorporated in the available data. 

The slight advantage of the first figure in (3.11a) in what regards the x^/d.o.f. makes one wonder 
that the difference is really caused by isospin breaking (in which case the value obtained from t decay should 
be rejected) or is due to random fluctuations of the data, as well as to the systematics of the experiments. 
The second explanation has in its favour that, if we include the spacelike data into the fit (but still use 
(3.9a)) the discrepancy is softened, and we get compatible results: 

to +/nQnA72N / 4 754 ±55 (e+e" + spacelike; X^/d-o-f. = 179/154) /qiik\ 
a(2.;i<0.8GeV) = |^g26±23 (r + spacelike; xVd-oi. = 112/93). (^-"^^ 

This last result allows us to draw the following conclusion: that part of the discrepancy between results 
obtained with e~^e~ and r decay is still of statistical origin, but also it would seem that part is genuine. 

In an attempt to take into account at least some of the isospin breaking effects, we have fitted 
simultaneously e+e~, r decay, both including spacelike data, allowing for different values of the mass and 
width of the rho (but keeping other parameters, in particular ci, C2, common for both e~^e~ and r fits). We, 
however, still use (3.9a). In this case we find convergence of the results; we have'' 



a(27r; t < 0.8 GeV^) = 4 779 ± 30, xVd-o.f. = 248/204; (e+e" + r + spacelike), (3.12) 

which is compatible (within errors) with both numbers in (3.11b). 

It is to be noted that, if we had not allowed for different masses and widths for the neutral and 
charged rho, we would have obtained, in this common fit, 

10^^ X a(27r; t < 0.8 GeV^) = 4 822 ± 30, xVd-o.f. = 264/206; {e+e~ + r + spacelike) 

i.e., a larger x^/d.o.f. and a value quite different from that obtained with only e+e" and spacelike data. 
So it would appear that allowing for different parameters for the neutral and charged rho really takes into 
account a good part of the isospin breaking effects. 

Finally, we take into account the kinematical effects of the tt^, tt" mass difference repeating the fit 
using (3.10b) now.^ The result of the fit with e+e~ data only is of course unchanged, but we reproduce it to 
facilitate the comparison and for ease of reference. We find what we consider our best results: 

10" X a(27r; t < 0.8 GeV^) = 4 774 ± 31, xVd-o.f. = 246/204; (e+e" + r + spacelike). 

10" X a(27r; t < 0.8 GeV^) = 4 754 ± 55, xVd.o.f. = 179/154; (e+e" + spacelike). ^^'^^^ 

We remark that the results for the evaluation including r decays are rather insensitive to the use of (3.10b), 
but what change there is, it goes in the right direction: the x^/d.o.f. has improved slightly, and the values 
for the anomaly including the r have become slightly more compatible with the figure obtained using e'^e~ 
data only. This makes us confident that most of the effects due to isospin breaking, both from u, d mass 
differences and from electromagnetic effects (about which we will say more in Subsects. 3.3 and 5.2) have 
already been incorporated in our calculation. The fit may be seen depicted in Fig. 2 for |i^7r|^, with timelike 
and spacelike data, and in Fig. 3 for the quantity vi in r decay. 

^ When evaluating a(27r; t < 0.8 GeV'^) we of course use the parameters nip, bo, 6i corresponding to p"; sec below. 

® For consistency we should also have taken the expresion k = ^{[t — (m^+ — m^o)^][t — (m^+ + m^o)^]y^^ , altered 
the threshold to t = (n^,r+ + ™-iro)^ for tau decay and allowed for different scattering lengths. We have checked 
that the effect of this on the contribution to a leaves it well inside our error bars; wo will discuss the results one 
gets in a separate paper. Note that it makes sense to still consider the same ci, C2 for e^e~ and tau decay as 
these parameters arc associated with G whose imaginary part vanishes below t = so ~ 1 GeV^ where the effects of 
isospin breaking should be negligible. 
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-0.2 0.2 0.4 0.6 0.8 1 1.2 

t (GeV^) 

Figure 2. Plot of the fit to \F^it)f, timelike (ref. 8) and spacelike (ref. 11) data. The theo- 
retical curve actually drawn is that obtained by fitting also r data, but the curve obtained fit- 
ting only e"'"e~ could not be distinguished from that drawn if we plotted it. A blowup of the fit 
in the spacelike region may be seen in Fig. 4. 



The /d.o.f. of the fits is slightly above unity; in next subsection wc will sec that including system- 
atic errors cures the problem. For example, just adding the systematic normalization error for the spacelike 
datal^^l gives a shift of the central value of 31 x 10~^^ and the x^/d.o.f. decreases to 152/153 for the evalu- 
ation with e+e" data only. The quality of the fit to the spacelike data is shown in Fig. 4, which is a blowup 
of the corresponding part of Fig. 2. 

The parameters of the fits are also compatible. We have, 

ci = 0.23 ± 0.02, C2 = -0.15 ± 0.03; bo = 1.062 ± 0.005, h = 0.25 ± 0.04 (e+e" + t + spacelike); 
ci = 0.19 ± 0.04, C2 = -0.15 ± 0.10; bo = 1.070 ± 0.006, bi = 0.28 ± 0.06 (e+e" + spacelike). 

(3.14) 

In the first line the parameters Ci, C2 are common for p^, p+. bo and bi vary very little; the ones shown 

correspond to the values of mpO, FpO as given below in Eq. (3.15). The values Bq = Constant, 6i = would 
correspond to a perfect Breit-Wigner shape for the p. Another fact to be mentioned is that including the 
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Figure 3. Plot of the fits to v\(t) (histograms), and data from t decay (black dots). 
Left: Aleph data. Right: Opal data. (ref. 10). Note that the theoretical values (his- 
tograms) are results of the same calculation, with the same parameters, so the differences 
between the two fits merely refiect the slight variations between the two experimental de- 
terminations. 




-0.05 

t (GeV^) 

|2 



-0.05 

t (GeV^) 



Figure 4. Plot of the fit to |Jir(t)| in the spacelike region. With only statistical errors 
(left) and including systematic experimental errors (right). 



corrected phase space factor (3.10b) helps a little to make compatible the parameters for both fits; if we had 

used (3.9a) wc would have obtained 

ci = 0.23 ± 0.01, C2 = -0.16 ± 0.03; ho = 1.060 ± 0.005, hi = 0.24 ± 0.04 (e+e" + r + spacelike). 
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An important feature of our fit is that the coefficients decrease with increasing order. This, together 
with the fact that the conformal variables w, z arc of modulus well below unity in the regions of interest 
{Ami <t<0.8 GeV^ for w, -0.25 GeV^ < t < 0.8 GeV^ for z): 

-0.57 < w; < 0.24, -0.38 < z < -0.02, 

ensures good convergence of the expansions. We have also checked that including extra terms in the expan- 
sions does not improve the quality of the fits significantly. 

Besides the results for the anomaly we obtain reliable determination of a set of parameters. We have 
those pertaining to the rho. 



rupo = 772.6 ± 0.5 MeV, FpO = 147.4 ± 0.8 MeV; 
mp+ = 773.8 ± 0.6 MeV, TpH- = 147.3 ± 0.9 MeV . 



(3.15) 



The figures arc in reasonable agreement with the Particle Data Group values''' given before. 

The value for the scattering length the fit returns is comfortably close to the one obtained from tttt 
phase shifts; we get 

al = (41 ±2) X 10^3 ^-3_ 

This value of al is slightly larger, but compatible with recent determinations based on an analysis of tttt 
scattering (ACGL) or chiral perturbation theory (CGL, ABT) that give (ref. 18) 

al = (37.9 ±0.5) x 10"^ m-^(CGL); al = (37 ±2) x 10"^ m"^ (ACGL); al = (38 ± 2) x 10"^ m"^ (ABT). 
Also from our fits we obtain the low energy coefficients of the pion form factor. 



t- 



(3.16) 



{rl) = 0.435 ± 0.002 fm^, c„ = 3.60 ± 0.03 GeV"^ (e+e" + r + spacelike); 

{rl) = 0.433 ± 0.002 fm^, c„ = 3.58 ± 0.04 GeV"^ (e+e" + spacelike). 

These figures are also compatible with, but much more precise than, the current estimates: I^^l 

(rl) = 0.431 ± 0.026 fm^ c„ = 3.2 ± 1.0 GeV"* . 

Another remark is that in all these fits we took to = 1.1 GeV^. The dependence of the results 
on this parameter, to, is very slight, provided we remain around this value. Thus, for example, if wc take 
^0 = 1.2 GeV^ the value of a(27r; t < 0.8 GeV^) only increases by 4 x 10^^^, and the global only varies 
by one unit. 

As further checks of the stability and reliability of our results wc mention the following two. First, 
we could, as discussed above, have imposed the more stringent values for aj as given in ref. 18. Now, if for 
example wc take, in accordance with (ACGL) in this reference, the value a} = (37 ±2) x 10~^ x m~^, instead 
of the value a\ = (38 ± 3) x 10^"^ x m^'' that follows from only experimental tttt data, the fit deteriorates. 
The fit returns the value a\ = (39 ± 1) x 10^'^ x m,^'^ for the scattering length, in (slightly) better agreement 
with the input; but we do not consider this an improvement as the global increases by 2 units. 

The corresponding value for the contribution to the anomaly changes very little, from the value 
(4774 ± 31) X 10^1^ (Eq. (3.13)) to (4768 ± 32) x 10"" now, i.e., a shift of only 6 x 10"" with a small 
increase of the error. Thus, the results are insensitive to a more stringent input for a} but, because the quality 
of the fit deteriorates, we still consider the result with the more relaxed input o} = (38 ± 3) x 10~^ x 
to be less biased. 

Secondly, we have not used the experimental phase shifts as input (except for the values of the 
scattering length). So, the values that follow from our expression (3.6), with the parameters given in (3.14), 
constitute really a prediction for $1 (t) . This can be compared with the existing experimental values for this 



It should be noted that the various determinations for nip reported by the PDG'^®' actually cluster around severaJ 
different values. 
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Figure 5. Our predicted phase 
TTTT shift, d\ (in radians), com- 
pared to the experimental values 
for the same (Solution 1 from Pro- 
topopescu et al., ref. 17). The ex- 
perimental errors are of the order 
of the size of the black dots. 



quantity, ["'^^1 a comparison that may be found in Fig. 5. The agreement is remarkable. The result one would 
have obtained if including the phase shifts in the fit may be found after Table 1. 

Before finishing this section we have to clarify the matter of the ui and ui — p contribution to 
a(27r; t < 0.8 GeV^). Our fits to e+e^ data have actually been made including in the function as given 
above, Eq. (3.2), a coefficient to take into account the w ^ 27r contribution. To be precise, we have used the 
expression 

1 + 0^ 2 • fr 

Ff{t) = F^{t) X m t im^i^ ^^ .^^^ 

1 + (J 

where the notation is obvious. We take from the FDCl^^l the values for the mass and width of the w, 

= 782.6 ± 0.1 MeV, F^ = 8.4 ± 0.9 MeV, 

and the fit gives a mixing parameter cr = (16 ± 1) x 10~^. 

As is known, this Goimnaris -Sakurai'^^^l parameterization is only valid for t ~ ra^ ^ and, in partic- 
ular, its extrapolation to f ~ is not acceptable. This effect is very small, less than one part in a thousand. 
However, to play it safe, we have also adopted the following alternate procedure: we have obtained a first 
approximation to F^^ by fitting the experimental data excluding the region 0.55 GeV^ < t < 0.65 GeV^. 
Then we have fitted only this region adding there also the uj piece, as in (3.17). The resulting value for 
a(27r; t < 0.8 GeV^) varies very little; it decreases by something between 2 and 12 x 10~^^, depending on 
the fit. We may consider this as part of the theoretical error of our calculation, to be discussed in next 
subsection. 

To finish this subsection, we present in Table 1 a comparison both with old results that also use 
analyticity properties, and a recent one (which does not).^° 

'^'^A different case is the analysis of A. Pich and J. Portoles, Phys. Rev. D63 093005, (2001), which also uses the 
Omnes equation method. This paper presents unfortunately a number of weak points. For example, the authors 
use an incorrect analyticity structure for 6l (their equations (8) and (A2), without left hand cut or inelasticity 
cut); they also employ a mere Breit Wigner to describe the phase in the rho region, while it is known that the rho 
deviates from this (e.g., our term bi). They also forget the cut at high energy in their equivalent of our G function. 
In what respects their results, the sitation is as follows. The value Pich and Portoles give (in units of 10~^^ and for 
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4 795 ±61 


Nl 


r + e+e" 


4 730 ± 100 


N2 


e+e" 


4 846 ±50 


CLY, AYW 


e+e" 


4 794 ±50 


CLY-IlI'*! 


e+e" + TTTT ph. shifts 


1771±:U 


T\'l 


r + c + v- 


4 754 ±55 


TY2 


only e+e~ 


Table 1 



Comparison of evaluations of lO'^^ x a{2Tv; t < 0.8 GeV^). 
Nl, N2 are in ref. 7. TY denotes our result here (statistical errors only for the e+e~ and spacelike data). 

The difference between the old CLY, AY and the new determinations is due to a large extent to 
the influence of the new Novosibirsk and NAT data which allow us in particular to obtain a robust result: 
the CLY evaluation used only 18 data points! In this respect, we note that, if we had included the tttt 
phase shifts in the fit (with also r decay data) we would have obtained 4781 ± 29 for this 27r contribution, 
i.e., a shift of only 7 units (as compared with a shift of 48 in the CLY-II evaluation). The value of the 
scattering length would be a} = (43 ± 3) x 10~^m~^ now. The corresponding x^/d.o.f., 276/227 with only 
statistical errors, is also good. This is an important proof of the stability of our results against introduction 
of extra data. (However, as explained above, we prefer the result without fitting phase shifts because of the 
model-dependence of the last). 

The difference between the results of Narison (N), who does not take into account the Fermi- Watson 
theorem or the spacelike data and TY, who do, is due in good part to, precisely, the influence of the spacelike 
data which also help reduce the errors. 

3.3. Systematic and theoretical errors for the pion form factor contribution 

Errors included in this work are divided into statistical and systematic. 

Evaluation of the statistical errors is standard: the fit procedure (using the program MINUIT) 

provides the full error (correlation) matrix at the minimum. This matrix is used when calculating the 
corresponding integral for a^, therefore incorporating automatically all the correlations among the various 
fit parameters. 

In addition, for every energy region, we have considered the errors that stem from experimental 
systematics, as well as those originating from deficiencies of the theoretical analysis. The experimental 
systematics covers the errors given by the individual experiments included in the fits. Also, when conflicting 
sets of data exist, the calculation has been repeated, and the given systematic error bar enlarged to encompass 
all the possibilities. 

In general, errors (considered as uncorrected) have been added in quadrature. The exceptions are 
explicitly discussed along the text. 

We next discuss the errors that stem from experimental systematics, as well as those originating from 
deflciencies of the theoretical analysis for the 27r contribution, in the low energy region 4m^ <t<0.8 GeV^. 
We start with the systematic errors of the data. They are evaluated by taking them into account in a new 

the contribution of 2n at energies below 1.2 GcV^) is 5110±60(PP), with a xVd-o.f.=33.8/21. This is substantially 
higher than all other results: 5027 ± 61 ( Narison, t and e+e~ data); 5044 ± 67 (Our result, only timelike r data), 
X^/d.o.f.=53/48, and 5004 ± 51 (Our best result, including e+e" and spacelike data), x^/d.o.f.=213/204. In the 
last three we have taken the piece 0.8 < t < 1.2 GoV'^, equal to 230 ± 5, from e+e~ data. The result of Davier 
and Hocker (ref. 5), using tau data only, is essentially like ours. No doubt the bias introduced by the wrong 
parametrization of the paper of Pich and Portoles is responsible for their discrepancy. 
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fit. In this way we find, in units of 10 and neglecting the mass differences corrections (i.e., using (3.9a) 
for tau data) 

Exp. Sys. = ± 40 {e+e- + r) 
Exp. Sys. =±66 (e+e"). 

To estimate the degree of correlation of the systematic errors pertaining to several experiments is a difficult 

task; wc choose to consider the full range from to 1. The error bars given cover all the possibilities. The 
X^/d.o.f. of the fit improves when taking these systematic errors into account to 

xVd-o.f. = 214/204 {e+er + r) 
xVd.o.f. = 145/154 (e+e"). 

The error given for the case in which we include the decays of the tau would be smaller, and the 
X^/d.o.f. would improve, if we used the correct kinematical formula for phase space, Eq. (3.10b); we would 
have obtained 

Exp. Sys. = ±30; xVd.o.f. = 213/204 (e+e" + r). 

In spite of this, we choose to accept the larger error (±40) as we feel that it includes residual effects of isospin 
breaking and electromagnetic corrections, different for the tau and e+e~ cases, that we will discuss at the 
end of this subsection. 

We discuss the systematic and theoretical errors in the higher energy regions in next section, but we 
mention here that the systematic error (4 x 10~^^) for 2-k between t = 0.8 and 1.2 GeV^ is added coherently 
to the lower energy 27r piece. 

In addition to this we have several theoretical sources of error. First, that originating in the ap- 
proximate character of the Gounnaris-Sakurai method for including the uj. This we estimate as discussed 
at the end of Subsect. 3.2, getting on the average ±7 x 10~^^. Then, the dependence of our results on to 
can be interpreted as a theoretical uncertainty, that we take equal to 4 x 10~^^. Composing these errors 
quadratically, we can complete (3.13) to 

n(9^- f < n S C.V^\ _ / 4 774 ± 31 (St.) ± 41 (Sys. +th.) = 4 774 ± 51 (e+e" + r + spacelike) 

^ ' - ~ \ 4754±55 (St.)±66 (Sys. +th.) = 4754±86 (e+e" + spacelike). ^ ' 

To finish this subsection we will discuss in some detail some matters concerning to radiative cor- 
rections and isospin breaking in as much as they affect the error analysis. We will start with the analysis 
based on e+e^ data. When evaluating the pion form factor we have used formulas, deduced in particular 
from unitarity and analyticity, that only hold if we neglect electromagnetic (e.m.) interactions. However, 
experimentalists measure the pion form factor in the real world, where the tt+tt" certainly interact electro- 
magnetically. Not only this, but the initial particles (e+e~) also interact between themselves, and with the 
pions. 

These last electromagnetic interactions, however, can be evaluated and they are indeed subtracted 
when presenting experimental data on F^^; the uncertainties this process generates arc estimated and in- 
cluded in the errors provided with the data. We will thus only discuss the uncertainties associated to e.m. 
interactions of the tt+tt" alone. These particles may exchange a photon, or radiate a soft photon that is not 
detected (see the corresponding figures in Subsect. 5.2). We may then define two quantities: Fi^\ which is 
the form factor wc would have if there were no cm. interactions; and which is the measured quantity, 

even after removal of radiative corrections for initial states or mixed ones. Actually, F^'^^'^^^ depends on the 
experimental setup through the cuts applied to ensure that no (hard) photon is emitted. 

Our formalism, as developed in Subsect. 3.1, applies to but we fit f},""^^\ Therefore, we are 

introducing an ambiguity 

r.(real) _ rn(0) 

which is of order a. 

The error induced by this ambiguity should be small. In fact, what enters into is the sum of the 
contribution of Fjr"^^^\ which is what we actually fit, and that of the state ■k~^tt~j, which can be obtained 
from the process 

e~'~e~ — > (7) — > 7r~''7r~7. 
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For this reason, we believe that the error due to the mismatch of Fi'^^'^^' and F^^^ is included in the errors to 
our fits hcrc;^^ the estimated error for the 7r"^7r~7 contribution, 9 x 10"^"'^, will be evaluated in Subscct. 5.2. 

Tau decay presents the same difficulties, and we expect a similar uncertainty as for e+e" collisions. 
But apart from the effect ^ discussed, it poses extra problems when relating it to -F^. To discuss 

this, we take first ^ rrid, a — 0; then we choose a ^ 0, but m„ = m^. Higher effects, proportional to 
a(m„ — rrici), o? and (m„ — m^'^ shall be neglected. 

For a = 0, the masses of 7r+ and 7r° become equal, but isospin invariance is still broken. This means 
that, in particular, the quantity 77"^ in (3.10a) is nonzero and therefore the experimentally measured v\ does 
not coincide with that in (3.10b). We expect this effect to be small, since it is of second order, (to^ — to„)^. 
If the scale is the QCD parameter A, then this will be of relative size 10~^; but other effects need not be so 
small. We have tried to take them into account by allowing for different parameters for p+, this produced 
a substantial shift, of about 40 x 10~^^ for a^. 

Then we set rn„ = and take e.m. interactions to be nonzero. Apart from the effects already 
discussed, this produces the mass difference between charged and neutral pions. This we took (partially) 
into account when using the modified phase space of (3.10b). The ensuing effect for turned out to be 
small, ~ 4 X IQ-^i. 

Remnants of a 7^ and m„ 7^ m,^ will likely still affect the determination of from e+e~ and 
r decay data; but we feel that accepting the systematic/theoretical error of 40 x 10~^^ covers the related 
uncertainty. 

4. Contributions to a(2)(h.v.p.) from t > 0.8 GeV^ The full a(2)(h.v.p.) 
4.1. The higher energy contributions, and the Stt contribution 

At higher energies we will get a substantial improvement over determinations based on old data^^"! because 
of the existence of very precise measurements from Novosibirsk 1^1 and Beijing, I^^l gathered in the last two- 
three years, which will help remove a large part of the existing errors. This is particularly true of the region 
up to f = 3 GeV^ which caused an important part of the total errors in pre-1998 calculations of a(^^(h.v.p.). 
We turn to it next. 

4.1.1. The region up tot = 3 GeV^ 

We consider first the contribution of two, three, four pion, . . . , and KK intermediate states for 
0.8 < t < 1.2 GeV^ In what follows n.w.a. will mean narrow width approximation, r.d.a. resonance 
dominance approximation (but not narrow approximation) and s.o.i.c. sum over individual channels. For 
the n.w.a. we use the standard formula. Denoting by FeeiV) to the width into e+e" of a vector resonance 
V with mass M, its contribution to a'-^^(h.v.p.) is given in this approximation by 

^ ^ ttM ^ ' 

The uncertainty on a(V) is calculated by Gaussian error propagation for the parameters in (4.1). In practice, 
it is dominated by the experimental error of the electronic width. 

Stt states, <t<\.2 GeV^ 

In the narrow width approximation one gets the cj, contributions; 

10" X a(37r;u;) =348 ± 13 

10^^ X a(37r;(/)) =62 ± 3, 

but this misses the region between a; and </>, and the interference effect just above the last 
experimental data. This gives 

10" X a(37r; t < 1.2 GeV^) = 438 ± 4 (Stat.) ± 11 (Sys.). 



(4.2) 



So we will use 



(4.3) 



In particular for the evaluation including r decay data, see below 
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FlGURE 6. Plot of the fit to the 

e^e~ —> Stt cross section up to 
t = 1.2 GeV^, with data from rcf. 9. 
Continuous Une: fit to CMD2 and SND 
data. Dashed fine: fit to CMD2 and ND. 



The fit to the Stt experimental cross section, with data from ref. 9, may be found in Fig. 6. The 
upper curve (continuous line in Fig. 6) is a fit to the CMD2 and SND data. We have used a Breit-Wigner 
parameterization for the lo and (f> resonances, plus a constant term and the exact threshold factor for Stt 
states. The x^/d-o.f. is 63/60; we consider this our central result here. The dashed curve fits instead the 
data from CMD2 and ND (Dolinsky et al, ref. 20); the quality of the fit is poorer (X^ /d.o.f. = 52/37). It 
fits better the region between the to and (j), but fails to reproduce the data beyond 1.06 GeV^. In fact, we 
include the second fit to estimate the corresponding systematic uncertainty; the small difference in terms of 
the integrals between the two fits, 8 x 10~^^, is included into the systematic error. 

2tt states, 0.8 GeV^ <t<1.2 GeV^ 



This 27r state contribution is 

10^^ X a(27r; t < 1.2 GeV^) = 230 ± 3 ± 4. (4.4) 

The evaluation of the contribution of the 27r state has greatly improved (with respect to older 
calculations) because of the information from recent Novosibirsk'^^ data on e+e^ 27r. We have fitted the 
experimental value of I^Vp with an expression l/{bt + a), a, b completely free parameters; the result of this 
fit may be seen depicted in Fig. 7. (A similar result is obtained if we extended our earlier calculation of 
Fjr{t) to t ^ 1.2 GcV^ by setting tg = 1.2; but we prefer the result based only on experimental data.) Of the 
two errors given for the 2tt contribution the first is statistical and the second, systematic, has been added 
coherently to the systematic error on the low energy 27r contribution, as discussed in Subsect. 3.3. 

The wiggle in Fig. 7 for t ^ m| is due to the interference of the decay (t> 27r. This is similar to 
the uj — p effect, and has been treated in a similar manner; we have incorporated it using the formulas and 
parameters given by Achasov et al.I^l The influence of this effect on is minute. 

KK states, 0.8 GeV^ < i < 1.2 GeV^ 



An important contribution is that of KK states. In the n.w.a., this is given by the </>: 

10^'^ X a{KK;(l)) = 2,2,2±9, (4.5) 
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FlGURE 7. Plot of the fit to 
\F^{t)f in the region 0.8 < t < 
1.2 GeV^ 



but this is a dangerous procedure here; the vicinity of the KK threshold distorts the shape of the resonance. 
We thus have to calculate this KK contribution directly from experiment. We have used two fitting pro- 
cedures. In the first, we fit simultaneously the K~^K~ and K^Kg data of Achasov et al.,[^l with the same 
parameters for the (!>. We get. 



10 



11 



a{K+K-; t < 1.2 GeV^) = 185.5 ± 1.5 (Stat.) ± 13 (Sys.) 



and 



10" X a{KLKs; t < 1.2 GeV^) = 129.5 ± 0.7. 

The quality of the fit, shown in Fig. 8, is good (X^ /d.o.f. = 84/82). 

In the second fitting procedure, we add the K^Ks data of Akhmetshin el al.,[^l obtaining the result 

10" X a{KLKs; t < 1.2 GeV^) = 128.4 ± 0.5 (Stat.) ± 2.6 (Sys.). 

The fit is now less good, but the integrals are essentially identical for both fits. Adding the KK results 
together we find 

10" X a{KK; t < 1.2 GeV^) = 314 ± 2 (Stat.) ± 13 (Sys.). (4.6) 

The systematic errors have been obtained repeating the fits, including now the systematic errors given by 
the experiments. 

We mention in passing that the ratio of contributions of K~^K~ and KlKs, 1.44, agrees well with 
the ratiol^^l 

r{4)^ K+K^) 
r{(i, ^ KlKs) 

Other statc^s: 4/1. 5" . //tt""" ■ ■ ■: 0.8 < f < 1.2 Gc-V^ 



1.46 ±0.03. 



The four pion contribution, including the quasi-two body state uiir, may be obtained from recent 
Novosibirsk data, or from the compilation of Dolinsky et al.^^l If we use the last we get 

10" X a(47r; t<l.2 GeV^) = 25 ± 4; 

if we fit the data of Akhmetshin et al.'^l we find 

10" X a(47r; t < 1.2 GeV^) = 18 ± 3. 
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Of the 477 contribution most is due to the wtt^ channel; only a small fraction (2.4 x 10 comes from the 
TT+TT" TT+TT" statcs. Wc take, for this 4tt contribution, the figure 

10^^ X a(4TT; t < 1.2 GeV^) = 20 ± 5, (4.7) 

which covers all possibilities. 

The five, six, . . . , pions as well as a; ^ + 2tt° contributions are very small. Altogether, they 

give 

10^^ X a(5TT, 67T, r/7T°7T°, • • • , t < 1.2 GeV^) = 4 ± 2. (4.8) 

We present the summary of our results in the important region 0.8 GeV^ < t < 1.2 plus the 3tt 
contribution below 1.2 GeV^ in the following Table 2: 



Channels 




Comments 


tt+tt" 


230 ± 3 ± 4 


0.8<t< 1.2 GeV^ 


37r 


438±4±11 


9ml < t < 1.2 GeV^ 


K+K- 


186 ± 2 ± 13 




KlKs 


128 ± 1 ± 2 




4-7T 


20 ±5 


including wtt" 


Multipion, ij + 27r, • • • 


4±2 





Total 



1 006 ± 19 



Syst. error for 27r not included 



Table 2 



Contribution to a^^^ of various channels up to t = 1.2 GeV^ (27r below 0.8 GeV^ not included). 
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Figure 9. Experimental data and vari- 
ous interpolations between the VMD cal- 
culation, for small t, and QCD for larger 
t. Prom AY, ref. 4; data from ref. 20. 



1.2 GeV^ <t<2 GeV^ 



We consider three determinations: 



270 ± 27 (Here) 

278 ±25 (s.o.i.c, CLYW) 

265 ± 22 (VMD+QCD; AY). 



(4.9) 



The first is obtained from a numerical integration of the data,!^"! with a parabolic fit. The method 
referred to as "VMD+QCD; AY", details of which can be found in the AY^ paper, consists in interpolating 
between a vector meson dominance (VMD) calculation for quasi-two body processes (wtt, pw, . . .), plus a 
Breit-Wigner expression for two-body channels (tttt, KK , . . . ) at the lower end, and perturbative QCD at 
the upper end, the interpolation being obtained by fitting experimental data (see Fig. 9). Because we want 
to present a result as model-independent as possible, we will take as our preferred figure that obtained here 
from experimental data: 



W X a(1.2 GeV^ <t<2 GeV^) = 270 ± 27. 



(4.10) 



2 GeV^ < i < 3 GeV^ 



240 ± 3 (yl) ± 3 (Cond.) (QCD) 

250 ±19 (N^l, r.d.a.; only e+e") 

r, 4.11 
276 ±36 (Nt^l, r.d.a.; e+e" + r) 

222 ± 5 (St.) ± 15 (Sys.) (J, exp. data). 

(J) here denotes an evaluation, integrating with the trapezoidal rule, of a compilation of experimental data 
supplied by F. Jegerlehner. 
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For the QCD calculations we take the following approximation: for n/ massless quark flavours, with 
charges Qf, we write 



a^Q^jl + ^ + (1.986- 0.115n/) (^)' 
+ \ - 6.64 - 1.20nf - O.OObni - 1.240 ~j{ i' ] (— 

I ' 3(x;/Q/)-i ^ TT ^ 

To this one adds mass and nonpcrturbative corrections. We take into account the 0{m?) effect for quarks 
with running mass rhi{t), which correct i?/"' by the amount^^ 

-3j2Ql^Ut) je + 28^ + (294.8 - 12.3n/) (^)'| f-i. 



Finally, for the condensates we add 



and 



2tt 



1 _ 23f^ 
27 ^. 



mi{ipiipi). 



We neglect the condensates corresponding to heavy quarks (c, b) and express those for u, d, s in terms of 
f^m^, fjc^K using the well-known PCAC relations. 

In the QCD calculation, the error labeled "Cond." is found by inserting the variation obtained 
setting quark and gluon condensates to zero, and that labeled A by varying the QCD parameter. For this 
parameter we take the recent determinations ^^^1 that correspond to the value 

a,,(M|) = 0.1172 ±0.003; 

to be precise, we have taken (in MeV, and to four loops), 

yl = 373±80, f < m^; /I = 283 ± 50, < i < m^; yl = 199 ± 30, i > m^. 

For the gluon condensate we take (agG^) = 0.07 GeV"^. 

The four evaluations give comparable results, with the r.d.a. ones larger, and presenting also larger 
errors. As proved by the reliability of QCD calculations of semileptonic t decays, a similar process in a 
similar energy range, we think perturbative QCD can be trusted here, so we select the corresponding value 
as our best result. We write thus 

10^^ X a(2 GeV^ < i < 3 GeV^) = 240 ± 6, (4.12) 

where we have added linearly the errors due to A and the condensates. 

As a verification of the reliability of the calculation, as well as the improvement it presents when 
compared with earlier determinations, in the rather involved energy range 0.8 < t < 3 GeV^ (including here 
the full 377 contribution), we compare our value here (adding, for the occasion, the channels w, ^ ^ 7r°7, r/7, 
see Subsect. 5.2) with that obtained by Narisonl''! who uses resonance dominance and s.o.i.c, and to the old 
CLYW evaluation, with s.o.i.c. and QCD: 



10" X a(0.8 GeV^ < t < 3 GeV^ +uj^3tt) =1 559 ± 34 (Here) 

10" X a(0.8 GeV^ < i < 3 GeV^ + w ^ 37r) =1 631 ± 46 (Narison, t + e+e") 



(4.13) 



10" X a(0.8 GeV^ < i < 3 GeV'^+to-^ 3tt) =1 675 ± 65 (Narison, e+e") 
10" X a(0.8 GeV^ < i < 3 GeV^ w ^ 37r) =1 618 ± 97 (CLY, e+e") 
The compatibility between the results, using different methods of evaluation for many pieces, is reasonable. 



^^The corrections of order may be found in the paper of Narison,'^' together with references. We have checked 
that the effect of this correction is smaller than the errors of the leading terms. 
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4.1.2. The region 3 GeV^ < t < 4.6^ GeV^ 

This is another region where the availabiUty of recent precise data'^^l in the neighbourhood of the cc 
threshold, previously poorly known, permits a reliable evaluation. As a byproduct, we get an experimental 
validation of QCD calculations. 

3 GeV^ <t<2^ GeV^ 



We use perturbative QCD here and get 

10^^ X a(3 - 2^ GeV^) = 120 ± 0.8 (yl) ± 0.8 (Cond.). 

22 GeV^ < f < 32 GeV^ 



We have now very good recent experimental data. So we present two determinations: 

IQi^ X a(2^ - 32 GcV^) = 200 ± 1 (A) (QCD) 

10" X a{2^ - 32 GcV^) = 210 ± 3(St.) ± 14(Sys.) (Exp. BES) 

We only give the error due to A here because that due to the condensates is negligible. When integrating the 
BES data we have used the trapezoidal rule. If instead we fitted a horizontal line, we would have obtained 

10^^ X a(22 - 32 GeV2) = 207 ± 2(St.) ± 13(Sys.) (Exp. BES) 

The BES[^2] purely experimental result and the QCD calculation are compatible, but one has to take into 

account the systematic errors of the first. This shows clearly the importance of systematic variations in e~^e~ 
annihilations data. We take as our preferred value for the sum of the two intervals that obtained from the 
QCD calculations: 

10" X a(3 GeV2 < i < 32 GeV2) = 320 ± 2 ± 1 = 320 ± 3. (4.14) 
32 GeV2 < t < 4.62 GeV2 



We give here the results in units of 10 We separate the contribution of the J/tjj, tjj', that we 
calculate in the n.w.a., and the rest. For the first we have, 

62.0 ±4.0 J/ij 
14.8 ±1.3 ip' 

For the remainder we have the following possibilities: 

91 ± 0.4 {A) uds; (QCD; 32 GeV2 < t < 4.6^ Ge\^) 
4.0 ± 0.4 tfj", t/>"', (N, r.d.a.) 

Total: 172 ± 4 (N; QCD±r.d.a.) 

91 ± 0.4 (yl) uds; (QCD; 32 GeV^ < t < 4.62 GeV^) 
46.8 - 28.6 ^ 38 ± 10 cc. (AY, NRQCD) 

(AY; QCD+NRQCD) 



Total: 206 ± 11 



54.7 ± 0.3 (yl) (QCD) : 3.02 < f < 3.72 
56 ± 0.3 ± 3 (Exp., BES) : 3.72 < ^ < ^ 

Total: 188 ± 4 ± 3 (Exp., BES+QCD) 
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Figure 10. Plot of BES experimental data and QCD for the u, d, s quarks (lower t) and the 
same plus NRQCD for the c quark contribution, from t = 3.74^ GeV^ to 4.6^ GeV^ at the 
right. Only systematic errors shown for experimental data. Statistical ones are even smaller. 



Here N refers to the paper of Narisonl^l, and AY is in ref. 4. BES are the experimental data from 
ref. 12. The first error for them is the statistical, the second the systematic one. 

This region merits a somewhat detailed discussion, as there is a certain controversy about it. We 
have made the calculation in three different manners. First, we separate the u, d, s quarks contribution, that 
can be evaluated using perturbative QCD. The contribution of the cc states is then evaluated saturating it 
by the ^jj resonances, in the r.d.a.; this is the result labeled (N, r.d.a.). This saturation procedure does not 
produce a good description. 

In a second method one separates also the u, d, s contribution; but the cc one is treated differently. 
If a resonance is below the channel for open charm production, which is set at t = 4m^ (with c the pole mass 
of the c quark), then it is treated as a bound state, in the n.w.a. Above cc threshold, one uses nonrelativistic 
QCD (see AY, ref. 4 and ref. 22). The two values reported above for such a calculation (AY, NRQCD) are 
for two values of the c quark mass: m,c = 1.750 GcV, in which case only the J/ip should be taken to be 
below threshold, and nic = 1.866 GeV and then both J/ip and ip' are to be added below threshold. This 
last gives the smallest number (28.6). The result of the calculation is taken as the average of both numbers, 
with half the difference as the estimated error. In Fig. 10 one can see the BESl^^l data and the predictions 
of QCD and NRQCD, the last for rric = 1.87 GeV. 

The third method, which is the one that yields our preferred number, 

10^^ X a(3^ GeV^ < t < 4:.&^ GeV^) = 188 ± 4 ± 3, (4.15) 

is obtained by using QCD for u,d,s quarks plus J/tp, tp' below t = 3.7^, and experimental data (BESl^^l) 
above that energy. 

All three methods give overlapping results, within errors, with the r.d.a. below experiment, and 
with an underestimated error, and with the NRQCD calculation reproducing better the data. This NRQCD 
calculation depends strongly on the mass of the c quark and, in fact, one can turn the argument backwards 
and predict rric by requiring equality with the experimental figure. If we do so, we find 

rUc ^ 1.89 GeV, 

a very reasonable estimate consistent with the two loop result, correct to 0{a^), which gives rric = 
1.866 ±0.20 GeV. 
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4.1.3. The region t > 4.6^ GeV^ 

The results from this region have not changed noticeably, but we give them for completeness. 
4.6^ GeV^ < t < 11.2^ GeV^ 



For the first T resonances, and in units of 10 , 

0.55 ±0.03 r 
0.18 ±0.01 T'. 

Then, 

88.8 ± 1.0 (A) udsc : (QCD), 4.6^ GeV^ < t < 11.2^ GeV^ 

0.22 ± 0.04 bb : (N, n.w.a.), T", T'", . . . 
0.53 ± 0.08 bb : (AY, NRQCD) 

Adding this, we get 

Total: 90 ± 1 (N; QCD+n.w.a.) 

Total: 90 ±1 (AY; QCD±NRQCD) 

The notation is hke for the c threshold region. The error in the (AY, NRQCD) evaluation is due to the 
error in the QCD parameter. A, and the b quark pole mass, for which we have taken^^l mb = 5.00±0.10 GeV. 
Both figures are essentially identical and we thus take 

10" X 0(4.6^ GeV^ < t < 11.2^ GeV^) = 90 ± 1. (4.16) 

11.22 Qgv2 <t^oo GeV^ 



The use of QCD is mandatory here. The contribution above tt threshold is negligible, so we calculate 
with n/ = 5 and get, 

10" X a(11.22 GeV^ < f ^ oo) = 21 ± 0.1 {A). (4.17) 

4.2. The whole a(2)(h.v.p.) 

Our final result for the 0(0^) hadronic contribution to is then 

10" X a(2)(h.v.p.) = il'f^mt + ^ + fr^-^ (4.18) 

^ \ 6 889 ±96 (e+e- + spacel.) ^ ^ 

To compare with other evaluations we have to add the contribution (43 ± 4) x 10^^^ of some of the 
radiative decays of the p, ui, (j) (see Sect. 5.2) that other authors include. This comparison is shown, for a 
few representative calculations,^^ in Table 3. 

^^A more complete list of evaluations, including some of the very earliest ones, may be found in the paper of Narison, 
ref. 7. 
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Authors 


10^^ X a(h.v.p.) 


Comments 


KNO 


7 nfiS + 1 74 


p~ nnlv 


CLY 


7 1 on + 1 1 fi 




CLY-II 




^ ~^ 1 c 1"! 1 ■77- TT r\Vi c Vi 1 c 


ADH 


7011 ± 94 




BW 


7 n9fi + 1 fin 




AY 


7 1 1 Q _(_ 1 no 




DH 


6 924 ±62 


e+e^ ± r 


J 


6 974 ± 105 




Nl 


7031 ±77 


e+e" + r 


N2 


7011 ± 117 


e+e" only 


TYl 


6 952 ±64 


e+e" ± r + spacel. 


TY2 


6 932 ±96 


e+e" ± spacel. only 



Table 3 



KNO: ref. 3; BW: ref. 3; J: ref. 6; CLY, CLY-II, AY: ref. 4; N: ref. 7; DH, ADH: ref. 5. 

If we had added also the other radiative contributions (tttt^, and the continuum hadron + 7, cf. 
Subsect. 5.2) we would have found the hadronic vacuum polarization piece, correct to order and a^, 

5. Higher order hadronic contributions 
5.1. Hadronic light-by-light contributions 

A contribution in a class by itself is the hadronic light by light one. So we split 

a(Other hadronic, 0{a^)) = a('One blob' hadronic, 0{a^)) + a(Hadronic light by light). (5.1) 

We will start by considering the last, given diagrammatically by the graph of Fig. 11. This can 
be evaluated only using models. One can make a chiral model calculation, in the Nambu-Jona-Lasinio 
version or the chiral perturbation theory variety, with a cut-off, or one can use a constituent quark model in 
which we replace the blob in Fig. 11 by a quark loop (Fig. 12). The result depends on the cut-off (for the 
chiral calculation) or on the constituent mass chosen for the quarks. After the correction of a sign error in 
the evaluation s of ref. 24 by M . Knecht and A. Nyffeler ( hep-ph/011105^ ), confirmed in M. Hayakawa and 
T. Kinoshita ( tiep-ph/0112102D we have 



10" X a(Hadronic light by hght) = 86 ± 25 Chiral calculation; BPP, HKS. (5.2a) 

Earlier calculations with the chiral model, using VMD to cure its divergence, gave (HKS, ref. 24) 

10" X a(Hadronic hght by light) = 52 ± 20 Chiral calculation (HKS). (5.2b) 

For the constituent quark model we use the results of Laporta and Remiddi.'^^^ The contribution to 
of light by light scattering, with a loop with a fermion of charge Qi, and mass larger than the muon 
mass, is 
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^4 fay 

where, to order {rrin/mi)'^, 



810 



9720 " 972000 



161 1„„2 16189 I 13 /-(o\ _ 161 „2 _ 831931 

48600 ^°6™ + 18''V<5J 

Taking constituent masses, 



m„,d = 0.33, rus = 0.50, rric = 1.87 GeV 



we find 



10^^ X a(Hadronic light by hght) = 46 ± 10 (Quark const, model) 
and the error is estimated by varying mu,d by 10%. 



(5.2c) 




Figure 12. The Ught by Ught hadronic 
correction in the constituent quark 
model. 
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One could also take the estimate of the 7r° pole from Hayakawa, Kinoshita and Sanda^^'^l and add 
the constituent quark loop, in which case we get 

lO" X a(Hadronic hglit by light) = 98 ± 22 (Quark const. mo(lel+ pion pole) : (5. 2d) 

one expects the chiral calculation to be valid for small values of the virtual photon momenta, and the 
constituent model to hold for large values of the same.^^ Indeed, almost a half of the contribution to 
a(Hadronic light by light) in the chiral calculation comes from a region of momenta above 0.5 GeV, where 
the chiral perturbation theory starts to fail, while for this range of energies, and at least for the imaginary 
part of (diagonal) light by light scattering, the quark model reproduces reasonably well the experimental 
data (see for example ref. 26 for a recent review of this). 

We will take here the figure, which comprises the relevant determinations, 

10" X a(Hadronic light by light) = 92 ± 20 (5.2e) 



5.2. Photon radiation corrections to the hadronic vacuum polarization 

The a('One blob' hadronic, O(ct')) corrections are obtained by attaching a photon or fermion loop to the 
various lines in Fig. 1. They can be further split into two pieces: the piece where both ends of the photon 
line are attached to the hadron blob, a(h.v.p., 7), shown in Fig. 13, and the rest. So we write, 

a('One blob' hadronic, 0{a^)) = a(h.v.p., 7) + a('One blob' hadronic, rest). (5.3) 

The last can be evaluated ^''1 in terms of the hadronic contributions to the photon vacuum polarization, 
finding 

10" X a('One blob' hadronic, rest) = -101 ± 6. (5.4) 

(Note, however, that this result has not, as far as we know, been checked by an independent calculation). 

The only contribution that requires further discussion is that depicted in Fig. 13, a(h.v.p., 7). In 
principle, this contribution can be evaluated straightforwardly by a generalization of the Brodsky-de Rafael 
method. We can write 

POO 

a(2)(h.v.p.) + a(h.v.p., 7) = / dtK{t)R^'^\t), (5.5) 

where 

n{2)i (T^^\e^e~ — > hadrons) + a^^^(e+e~ — > hadrons) + u^^\e^e~ — > hadrons; 7) 
^^'^ a(0)(e+e- ^/x+M") ' 

The notation means that we evaluate the hadron annihilation cross section to second order in a, and we add 
to it the first order annihilation into hadrons plus a photon. 



'^''Strictly speaking, one would also need large momentum of the external photon to get a really trustworthy evaluation 
with the constituent model. 
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/ . 7. hard 

y, soft 



Figure 14A. Diagrams included in the pion Figure 14B. Diagrams not 

form factor. included in the pion form factor. 



For energy (t) large enough, this can be calculated with the parton model, and leads to a correction 
f (E/ Q)/Y.f Q/)"/?!" times the parton model evaluation. Taking then t > 1.2 GeV^, this is (0.76 ± 0.04) x 
10~^^. The error is that due to A and the masses of c, h quarks. We have excluded the contribution of the 
radiative decays of the J/ f/", ij:'T, T' resonances since we have taken these into account already (we took the 
full e^e~ width for them). 

Then comes the contribution of small momenta, t < 1.2 GeV^. We start by discussing the process 
involving two pions. In our determination in Sects. 3, 4 of a(^^(h.v.p.), we made calculations by fitting the 
experimental cross section e+e~ — > tt+tt", which specifically excludes radiation of hard photons (hard pho- 
tons defined as those that are identified experimentally). Diagrammatically, this means that our evaluations 
of Sects. 3, 4 included the diagrams of Fig. 14 A (where a soft photon is one that is not detected), but not 
those of Fig. 14B (radiation of a hard photon). So, we have to include this radiation into a(h.v.p., 7). This 
can be easily done if we consider this region to be dominated by the rho, hence we approximate 

a^^\e'^e~ — > hadrons; 7) ~ ct'-^-' (e"'"e~ — > (p) — > 7r"'"7r~7). 

The last can be evaluated in terms of the branching ratio for the decay p — > 7r+7r~7, which is indeed measured 
experimentally (see the review of Dolinsky et al., ref. 20) from the reaction e+e" — > p — > 7r"'"7r~7. In the 
narrow width approximation for the rho, the contribution to is 

rjp ^ 77+77-7) 3ree(p)i^(m^) 

In this way, we find 

10^^ X a(h.v.p., 77+77-7) = 45 ± 7 (n.w.a.), (5.6b) 

and the error is that induced by the experimental error in the width r{p 77+77-7). 

We will elaborate a bit more on this contribution. The final state interaction of the 77+77- in the 
state 77+77-7 is very strong. The pions are produced in an S-wave, which presents a wide enhancement I^^l in 
the energy region E„+„- ~ 0.6 ± 0.2 GeV. However, this is only a small part of the contribution to the rate 



- 26 - 



-PRECISION DETERMINATION OF THE PION FORM FACTOR AND CALCULATION OF THE MUON Q — 2- 



for ■KIT'). According to Dolinsky et al.,P°l pp. 126 ff., most of the effect would be due to Bremsstrahlung by 

the pions. Above procedure to estimate this, in terms of the 7r+7r~7 decay of the p would be exact only if 
the experimental cuts made for identifying this decay, and to measure the pion form factor were the same. 
A more accurate procedure is as follows. We write 

.1.2 

o(7r+7r"7, t < 1.2) = / At K{t)R^+^-^{t) (5.7a) 

where 

R^^^-^{t) = B{t,E^)R^^^-{t) 
and the Bremsstrahlung factor B is given by'^^l 



m = km ( - k] log i±| - \km{\t''' -k)-k 



2iV2 ; i-^, 



Here ^ = [{km — k)/{\t^/'^ — k)]^/"^, km = (t — 4m^)/2t^/^ is the maximum energy of the photon and, finally, 
is the minimum energy for photon detection. 

To evaluate E^, we have to look at the setup of experiments measuring the pion form factor. 
Typically, one takes that no (hard) photon has been emitted when the angle between the pion momenta 
differs from tt by less than a small given amount, r]a. The energy cut is, in this case, 



E-r ^ Vo- 

The effective r]o depends on the specific cuts made in each experiment; those in ref. 8 are covered if we take 
?7o = 0.15 ± 0.05. Using this we find the result, for this range of rio, of 

10^^ X 0(77+77-7, t < 1.2) = 46 ± 0.5 ± 9. (5.8) 

The first error corresponds to the error in the integral of l-FVP and the second is induced by the 
dispersion in the values of t]o of the various experiments. (5.8) is practically identical to the n.w.a. result, 
Eq. (5.6b) (which of course is a satisfactory result). The reason for this agreement is that, when detecting 
7r+7r^7, the energy cut made, E~^ = 50 MeV (Dolinsky et al., ref. 20), turns out to be very similar to the 
average energy cut made when measuring the pion form factor. 

A similar analysis ought to be made, in principle, for other radiative intermediate states like 877 + 7 
and KK + 7, which can be estimated in terms of the corresponding decays of the w and 4>, but they give a 
contribution below the 10"^^ level and we neglect them. 

The lowest energy contributions to cr^^^e'^e~ — > hadrons;7) are those of the intermediate states 
7r°7 and 777, Fig. 15. At energies below the rho mass, one can evaluate the first (the only one that gives 
a sizable contribution) by relating the process to the decay 7r° — » 27. We write an effective interaction, 
corresponding to the vertex factor in the Feynman rules of 



k' 

with it 



2^ C/HTO/Jf^l f^2 i 



r(7r" ^ 27) = 
^ 25677 

so that = 4.6 x 10"^. Then, with e the electron charge. 
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tion to a(h.v.p., 7). 



This gives a very small contribution to a^, about 0.76 x 10~^^ if we integrate up to f^/^ ~ 0.7 GeV, and 
0.96 X 10~^^ if we go to i^/^ ~ 0.84 GeV (the integral only grows logarithmically). We only integrate up to 
the rho, i.e., to t^^"^ = 0.7 GeV with this pointlike model. 

Around the p region we have to take into account the excitation of this resonance, which produces 
the corresponding enhancement. This piece can be obtained in terms of the radiative width p 7r°7. More 
important is the w 7r°7 process which gives (33 ± 2) x 10~^^. Likewise, the contribution of the 777 state 
is evaluated in terms of the decay 777. Finally, the contribution from tt^tt^j is taken from ref. 29. 
Collecting all of this, we get 

10" X a(h.v.p., p 7r°7) = 4 ± 1 

10" X a(h.v.p., p^rj'j)^ 1.1 ± 0.4 
10" X a(h.v.p., w^7r°7) = 33±2 (5.9a) 

10" X a(h.v.p., (/) ^ 777) = 5 ± 1; 

Total : 43 ± 4 

we have included the lower energy contribution of ir^'y into a(h.v.p., p 77*^7) and, because we are relying 
on models, we added the errors linearly. For the 77777 states. 



10" X a(h.v.p., 7r+7r"7) = 46 ± 9 
10" X a(h.v.p., 7r°7r°7) = 2 ± 0.3; 

and, for the high energy piece, 

10" X a(hadrons + 7, t > 1.2 GeV^) =1 ± 0.5. 



(5.9b) 



(5.9c) 



Adding other contributions that are below the 10 level (e(700)7 ~ 0.7 x 10 etc.) we get the total 
effect of the states hadrons + 7, 

10" x a(hadrons + 7) = 93 ± 11. (5.10) 
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Channel 


Energy range 


Method of calculation 


i;rini"T'i rill "t" inn Vr\ 

10" X o(h.v.p.) 


TT+TT" 


t < 0.8 GeV^ 


Fit to e+e" + t + spacel. data 


4 774 ±51 




0.8 < t < 1.2 GcV^ 


Fit to cxp. e^e~ data 


230 ±5 


Stt 


t < 1.2 GeV^ 


B.-W. + const, fit to e^e~ data 


438 ± 12 


2K 


t < 1.2 GeV^ 


B.-W. + const, fit to e+e" data 


314 ± 13 


4tt, 5tt, tjtt, . . . 


t < 1.2 GeV^ 


Fit to e^e~ data 


24 ±5 


Inclusive 


1.2 GeV^ <t<2 GeV^ 


Fit to e^e~ data 


270 ± 27 


j/V^, V''; r, r 




N.w.a. 


77.5 ± 4.4 


Inclusive 


3.7^ GeV^ < t < 4.6^ GeV^ 


Fit to experimental e^e~ data 


56 ±3 


Inclusive: i/c/.s 


2 Gv\- < f < :>>-- G('\'- 


P(-rtuii)ativ(^ OCD 


615 ± 9 


Inclusive; udsc 


4.6^ GeV^ < t < 11.2 GeV^ 


Perturbative QCD 


89 ±1 


h quark 


10.1^ GeV^ < t < 11.2^ GeV^ 


Nonrelativistic QCD 


0.5 ±0.1 


Inclusive 


11.2^ GeV^ < i < oo 


Perturbative QCD 


21 ±0.1 


7+ hadrons 


Full range 


Various methods 


93 ±11 



Table 4 



Summary of contributions to a(^+^\ with what we consider the more reliable methods, as used in the present 
work. "B.-W. + const." means a Brcit Wigner fit, including the correct phase space factors, plus a constant; 
note that only for the four narrow rc;sonanccs J/V", V''; ^' use the n.w.a. The errors arc imcorrelated 
except those for QCD calculations (that have to be added linearly) and those for the 2?: states, for whose 
treatment we refer to the text. The errors given include statistical, systematic and (estimated) theoretical 
errors. For the details of the final states 7± hadrons we refer to Eqs. (5.9). 

6. Conclusions 

We present first, for ease of reference, Table 4 with a summary of the results obtained for a(h.v.p.) in the 
previous sections. 

Taking into account all contributions, and errors, we complete the best values for the h.v.p. piece, 
and the whole hadronic part of the anomaly: 

10" X a(2+3) (h.v.p.) = I ! 2S i 1 ^ ^ Tr^'^ (6-2) 

^ ^' \ 6 982 ±97 (e+e~ ± spacel.), ^ ' 

and adding the other radiative and light-by-light corrections, 

inii /-XT J • ^ [6 993 ±69 (e+e~ ± r ± spacel.) ,„ 
10 xa(Hadromc) = |gg^3^gg (e+e" ± spacel.) ^^'^^ 

Eq. (6.3) is of course the main outcome of the present paper. Because, even after adding systematic and 
theoretical errors the evaluation including r decay data is more precise, we may consider it to provide the 
best result for a (Hadronic) available at present. 
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Figure 16. Theoretical results on o(Hadronic) x 10 and experiment. 



J: ref. 6: Nl: ref. 7, data from e^e + t. N2: id., e^e only. Tl, T2: this paper with data from e^e + t or data 
from e^e~ only, respectively (including syst. and th. errors). 

We can add to (6.3) the pure electroweak contributions and present the result as the standard model 
prediction for afj,: 

10^^ X = 116 591 849 ± 69 (e+e" + r + spacel.) (6.4) 

We will next compare our results with other recent evaluations of the same quantities in Fig. 16, 
together with the experimental band. (They are shown incorporating the contribution of the 7r+7r~7 and 
TT*^ 77*^7 channels). 

A further point to emphasize is the importance of using, in the low energy region, parametrizations 
of F^(f) compatible with unitarity and analyticity. Only in this way we can incorporate data on F^(f) for 
spacelikc t into the fits. As discussed in Subscct. 3.2, these data not only provide a substantial shift for 
(of 39 X 10^^^ in the evaluation with e+e^ data only) but, by so doing, allow compatibility of these with 
the results from r decay, hence allowing a combination of the two in a meaningful way: this permits an 
important reduction of the errors of the calculation.-"^^ 

To finish this section, we can add a few words on prospects for improvements. In our view they 
are rather dim in the sense that it is not easy to see how one could get an error estimate clearly below the 
70 X 10~^^ mark, when taking into account systematic and theoretical errors. In fact, the central values 



^*Or, put conversely, not using data at spacelike t for implies a hidden error of about 40 x 10 
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have moved little, and the errors have not improved much, since 1985. It is true that experiments planned 
or in progress can clear further the region between 1.2 and 3 GeV^. However, a serious improvement of 
the very important low energy region for tttt is unlikely: as our evaluations show, one can get a fit to 
all data relevant for the hadronic component of o^, with a x^/d-o.f. of unity and verifying all theoretical 
constraints, with an error of at least 51 x 10""'^^, already for t < 0.8 GeV^ (sec Eq. (3.19)). In this respect 
the improvement obtained by adding r decay data, although not negligible, is minor: statistical errors are 
smaller, but theoretical ones are increased. 
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